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BASIC PROPERTIES OF NONLINEAR STOCHASTIC 
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q ! By Carlos M. Mora 1 and Rolando Rebolledo 2 

Universidad de Conception and Pontificia Universidad Catolica de Chile 

Mh| The paper is devoted to the study of nonlinear stochastic Schro- 

■^E . dinger equations driven by standard cylindrical Brownian motions 

(NSSEs) arising from the unraveling of quantum master equations. 
Under the Born-Markov approximations, this class of stochastic evo- 
lutions equations on Hilbert spaces provides characterizations of both 
f*S, ' continuous quantum measurement processes and the evolution of 

H , , quantum systems. First, we deal with the existence and uniqueness of 

regular solutions to NSSEs. Second, we provide two general criteria 
for the existence of regular invariant measures for NSSEs. We apply 
C3 , our results to a forced and damped quantum oscillator. 

1. Introduction. 

£> ■ 1.1. Objectives. In this work we focus on the nonlinear stochastic 

Schrodinger equations driven by standard cylindrical Brownian motions that 
describe open quantum systems under the Born-Markov approximation. 
More precisely, this paper concentrates on stochastic evolution equations 
on a separable complex Hilbert space ([),(•,•)) of the form 
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2S ■ (1) X t = X + V) / (L k X s - Re(X s ,L k X s )X s ) dW* 



V i + / (GX S + ]T(Re(X„ L k X s )L k X s - ±Re 2 (X s ,L k X s )X s ) da, 
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2 C. M. MORA AND R. REBOLLEDO 

where ||Xo|| = 1 a.s., W 1 , W 2 , . . . , are real valued independent Wiener pro- 
cesses on a filtered complete probability space (£l,$,($ t )t>o^) and G,L\, 
L 2 , • • • , are linear operators in f) with Dom(G) C Dom(Lfc), for any k G N, 
such that 

oo 

(2) 2Rc(x,Gx) + ^||L fc x|| 2 = 

fc=i 

whenever x € Dom(G). We are mainly interested in two problems: (i) ex- 
istence and uniqueness of the solution of (1); and (ii) existence of regular 
invariant measures for (1). 

1.2. Motivation. The primary motivation for this article is to develop 
the mathematical modeling of infinite-dimensional open quantum systems 
by means of stochastic processes. Our interest in the study of the nonlinear 
stochastic Schrodinger equation (1) arises mainly from the following three 
applications. 

First, (1) provides characterizations of the evolution of density operators 
and quantum observables. The classical model of open quantum systems 
consists of a small quantum system, whose state space is f) and its internal 
dynamics is governed by the Hamiltonian H, interacting with a heat bath 
or reservoir. Under the Born-Markov approximation, the time evolution of 
the density operators (positive operators in [) with unit trace [16, 56]) obeys 
the following quantum master equation in Lindblad form: 



Pi 



Po + J [g Ps + Ps G* + Y, L kPsL* k J da. 



Here Li,L 2 , ■ ■ ■ , take into account the effect of the environment and G is 
the effective Hamiltonian, that is, G = —iH — ^ J2k^=i -^it-^fc ( see [13, 28] for 
more details). Let us present a concrete physical model. 

Example 1 (Forced and damped quantum oscillator). Let t) = / 2 (Z + ). 
Suppose that (e n ) ng z+ is the canonical orthonormal basis on the space 
P(Z + ). The linear operators a) and a are defined in {x £ l 2 CL + ) 'J2n>o n \ x n\ 2 < 
oo} by a' e n = \Jn + le n+ i and 



ae r . 



(0, ifn = 0, 

Iv^n-i, ifn>0. 



Set N = a^a, the number operator. 

Choose H = i[ii(a) - a) + f3 2 N + f3 3 (a^) 2 a 2 . The interaction of the small 
system with the reservoir is simulated by L\ = a\a, L 2 = a 2 a) , L% = a^N, 
L4 = 04a 2 , L5 = «5(a^) 2 , and Lq = cxqN 2 , where Pi,f3 2 ,P3 are real numbers 
and «i, . . . , «6 are complex numbers. Consider L^ = for all k>7. 
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In Example 1, () represents the state space of a single mode of a quantized 
electromagnetic field and the vectors e n , with n G Z+, provide the energy 
levels of the system. Because a destroys a photon, Li^L^ model photon 
emissions. The operator a) describes the creation of a photon (see, e.g., 
[16, 56]). 

Using (1), we can obtain a probabilistic representation of pt- Indeed, it is 
to be expected that pt = E,\X t )(X t \ (see, e.g., [6, 31, 51]). In Dirac notation, 
|x)(x|, with x G rj, stands for the linear operator defined by |x)(x|(z) = 
{x,z)x for any zgf). Therefore, the probability that a measurement finds 
the system in the pure state x at time t > is E|(x,Xi)| 2 assuming that 
E|Xo)(A^o| is the initial density operator and x is a vector of () of norm 1. 
On the other hand, the value of the observable A at time t %(A) satisfies 
E(X t ,AXt) =E(X ,T t (A)X ) (see, e.g., [3, 34, 35]). Recall that quantum 
observables are represented by self-adjoint operators in rj. 

Second, X t is interpreted as the evolution of the state of a quantum system 
conditioned on continuous measurement (see, e.g., [4, 8, 19, 30, 51, 59]). For 
instance, the following example describes the simultaneous monitoring of 
position and momentum of a quantum system whose evolution is governed 
by the Hamiltonian H (see, e.g., [32, 57]). 

Example 2 {Continuous quantum measurement process). Let f) = L 2 (M, C). 
The position operator Q: rj — ► rj is given by Qf(u) = uf(u) for any mGK. 
The momentum operator P : rj — > rj is defined by P = —iD, where Df stands 
for the weak or distributional derivative of a function / G ff 1 (]R,C). Then, 
in the setting of (1) choose L\ = ^Q and L2 = naP, with k, o G ]0, +00 [. Set 
L k = for all k > 3. 

Third, (1) plays a relevant role in the numerical simulation of the time 
evolution of quantum systems (see, e.g., [31, 46, 51, 58]). In fact, using (1), we 
can overcome the difficulties arising in the direct numerical integration of the 
master Markov equations (i.e., quantum master equations and Heisenberg 
equations of motion) in Lindblad form when the dimension of the Hilbert 
space required for numerical computations is large (see, e.g., [13, 17, 51, 58]). 

Another motivation for this paper came from investigations in which (1) 
represents objective (independent of any observer) trajectories for quantum 
systems (see, e.g., [31, 51] and the comments of [61]). Furthermore, (1) ap- 
pears in the quantum filtering theory (see, e.g., [8, 10, 11]). This interesting 
application has been developed by Belavkin in the framework of quantum 
stochastic calculus [23, 33, 43, 49]. 

1.3. Outline of the paper. If rj is finite-dimensional and at most a finite 
number of L k are different from 0, then the existence and uniqueness of a 
strong solution of (1) can be obtained by means of customary techniques 
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employed in stochastic differential equations with locally Lipschitz coeffi- 
cients (see, e.g., Lemma 5 of [46]). In [5], Barchielli and Holevo covered the 
existence of a weak solution of (1) when G and L±,L2, ■ ■ ■ , are bounded 
operators. 

Using the linear stochastic Schrodinger equation 

ft °° ft 

(3) ¥*(£)=£+ / G^COete + E / !*¥>-(£) <W*. 



ft °° ft 

¥*(£)=£+/ GMOds + Y, U 
Jo k=1 Jo 



Girsanov's theorem and Ito's formula, Gatarek and Gisin treated in [29] the 
existence of solutions of (1) in two examples: (a) H = 0, L\ self-adjoint and 
L k = for any k > 2; (b) f) = L 2 (R, C), H = P 2 , L\ = Q, with P, Q defined 
as in Example 2, and L^ = whenever k>2. Exploiting Li is a self-adjoint 
operator, [29] also verified the pathwise uniqueness of solutions of (1) in 
(a) and (b). From the work of Kolokoltsov (see, e.g., [41]) it may deduce 
a generalization of the results of [29] for some multidimensional versions of 
the case (b). In [34, 35], Holevo sketched the proof of the existence of a 
weak solution of (1) in situations where, for any t > 0, (||¥> s (£)|| ) a e[o,t] 1S 
uniform integrable and E||93t(^)|| 2 =E||£|| 2 . To the best of our knowledge, 
the question of uniqueness of solutions of (1) is still unanswered (even in the 
case where G, L\, Li-, ■ ■ ■ , are bounded operators). 

In order to provide a sound basis for the study of infinite-dimensional 
open quantum systems by means of (1), Section 2 establishes the existence 
and uniqueness in distribution of the regular solution of (1) under a nonex- 
plosion condition on G and L\, Lz, . . . , similar to that used by Chebotarev 
and Fagnola to prove the Markov property of minimal quantum dynamical 
semigroups in [14] (see also [15, 23]). To this end, we combine the method 
of drift transformation (see, e.g., [36, 44, 53]) with the subtleties of (1). 

In recent years the large time behavior of quantum Markov semigroups has 
been the subject of a number of investigations (see, e.g., [2, 24, 25, 26, 60]). 
However, general results on the large time behavior of stochastic Schrodinger 
equations on infinite-dimensional Hilbert spaces of the types (1) and (3) 
are lacking in the literature. It is worth pointing out that Kolokoltsov [40] 
obtained that Xt is asymptotically similar as t — ► oo to a Gaussian function 
with certain time-dependent random parameters in case f) = L 2 (M,C), H = 
HP 2 /2, L\ = Q and L% = L3 = ■ ■ ■ = 0, where h denotes the Planck constant 
and P, Q are as in Example 2. For this purpose [40] uses an explicit solution 
of (3) (see, e.g., [9]). 

Since \\Xt\\ = 1 a.s., standard techniques of finite-dimensional stochastic 
processes yield the existence of stationary distributions for (1) when the 
dimension of f) is finite (see, e.g., [22]). A few attempts have been made to 
develop sufficient conditions for the uniqueness of an invariant measure for 
(1) in case dimrj < 00 (see, e.g., [7, 39]). 
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Section 3 deals with the basic problem (ii) . Under underlying assumptions 
similar to those in [24] , Section 3 states the existence of an invariant proba- 
bility measure T for (1) satisfying L \\Az\\ 2 T(dz) < oo, with A belonging to 
certain class of linear operators. This regularity property leads to the exis- 
tence of a quantum stationary state p^ for which the trace of Ap^ is well 
defined for a broad class of unbounded operators A. Here, our main criterion 
for the existence of regular invariant measures is based on characteristics of 
the operators G and L^. Moreover, this criterion involves the existence of a 
Lyapunov function inherent in the open quantum systems set-up. As a by- 
product, Section 3 provides the continuous dependence of the distribution 
of X on the initial data and the Markov property of X. 

Finally, Section 4 illustrates our main results with a forced and damped 
quantum oscillator. We select Example 1 as a model problem due to the 
role played by the one-dimensional quantum harmonic oscillators in the 
understanding of quantum systems (see, e.g., [16, 56]). 

2. Existence and uniqueness. This section provides a detailed study of 
the problem (i) in a context similar to that of [14] (see also [15, 23, 35]). We 
begin by specifying notation. 

2.1. Notation. Throughout this paper, the scalar product (•, •) is linear in 
the second variable and anti-linear in the first one. Furthermore, we assume 
that (e n ) ng z + is an orthonormal basis of [). Let rj n be the linear manifold 
spanned by eo, • • • , e n . We define P n : t) — > f) n to be the orthogonal projection 
of f) over \] n . 

Let A be a linear operator in f). Then Dom(A) stands for the domain of 
A and A* denotes the adjoint of A. The function 7ryi : f) — > fj is defined by 

( \_/ x > if x £ Dom(A), 
^W-jo, if x^Dom(A). 

Suppose that C is a self-adjoint positive operator in f). For any x,y £ 
Dom(C), we set (x,y)c = {x,y) + (Cx,Cy) and ||x||c = \/(x,x)c- Since C 
is a closed operator, (Dom(C), (•, -)c) is a Hilbert space. We write A instead 
of A o ttq as soon as the context avoid any confusion. 

The Borel u-algebra of the topological space E is written *B(E). 

2.2. Main result. In this section we suppose the existence of a reference 
operator C with the following properties. 

Hypothesis 1. The linear operator C : fj — >- fj is a self-adjoint positive 
operator such that: 

(Hl.l) Dom(C)cDom(G)nDom(G*). 
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(H1.2) For any n G Z+, t) n C Dom(C) and ££Li ||L£e n || 2 < oo. 
(HI. 3) There exist constants a,/3 G [0, +oo[ satisfying 

oo 

2Re{Cx,CP n Gx) + J^ \\CP n L k x\\ 2 < a(||Ca:|| 2 + \\x\\ 2 + (5) 
fe=i 

for any n G Z + and x G f) n . 
(H1.4) For all s G Dom{C), sup neZ+ ||CP n x|| < ||Cx||. 

Remark 1. Hypothesis 1 is a nonexplosion condition inherent in the 
open quantum systems context (see, e.g., [15, 23]). It applies to a broad 
range of applications as, for instance, models for heavy ion collisions [14], 
quantum oscillators (see, e.g., Section 4.3 of [23]) and quantum exclusion 
processes [47]. 

The notion of a regular solution of (1) given below is deeply inspired on 
the concept of a smooth classical solution of a partial differential equation. 
Loosely speaking, Definition 1 replaces the partial derivatives of a function 
by operators C satisfying Hypothesis 1 in order to describe the smoothness 
of a solution of (3). 

Definition 1. Let C satisfy Hypothesis 1. Suppose that T is either 
[0,+oo[ or [0,T], provided Te [0,+oo[. We say that (0,3, (3 r t ) teT ,P, 
(Xt)t£j,(Wt) t £j) is a solution of class C of (1) with initial distribution 
6 on the interval T if: 

• W 1 , W 2 , . . . , are real valued independent Brownian motions on the filtered 
complete probability space (Q,^,(^t)teT^)- 

• (X t )teT is an rj-valued process with continuous sample paths such that 
the law of Xq coincides with 9 and P(||Xt|| = 1 for all t G T) = 1. 

• For every t G T, Xt G Dom(C) P-a.s. and sup se r 0it i Ep||CX s || 2 < oo. 

• P-a.s., for all £ G T, 



ft °° ft 

(4) X t = X + G(X s )ds + y2 / L k (X s )dW : 

Jo f-r, Jo 



k 

s ' 

fc=r' n 



where, for any y G f), 

G(y)=Go7T C (y) 

OO 

+ ^(Ke{y,L k oTT C (y))L k oTT C {y) - lRe 2 (y,L k oTT C (y))y) 
fc=i 

and for any y G f) and k G N, 
(6) L k (y) =L k oTT C (y) - Re(y,L k o n c (y)}y. 
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For abbreviation, we simply say (P, (X t )teT, (Wt)tej) is a C-solution of (1) 
when no confusion can arise. 

The following theorem asserts the existence and uniqueness of the weak 
(in the probabilistic sense) regular solution of (1). 

Theorem 1. Let C satisfy Hypothesis 1. Suppose that 9 is a proba- 
bility measure on 23(b) such that 9(Dom(C) n {x G h : ||x|| = 1}) = 1 and 
L \\Cx\\ 2 9(dx) < oo. Then, (1) has a unique C-solution (Q, (Xj)t>o, (Bt)t>o) 
with initial distribution 8. Furthermore, for any t>0, we have 

(7) E Q ||CX t || 2 < exp(at)(EQ||c7X || 2 + ta(E Q ||Xo|| 2 + /?)). 

Theorem 1 shows that Example 2 equipped with a standard Hamiltonian 
is mathematically sound. 

Corollary 1. Let assumptions of Example 2 hold. Select H = aP 2 + 
(3Q 2 , with a > and (3 G ]R, and consider the self-adjoint operator [in L 2 (M, C)] 
N = {Q 2 + P 2 — -0/2. Suppose that 6 is a probability measure on QS(L 2 (M, C)) 
such that, for a given p G N, 9(Dom(N p ) D{x Gi 2 (K,C): \\x\\ = 1}) = 1 and 
f. \\N p x\\ 2 8(dx) < oo. Then, (1) has a unique solution of class N p with initial 
distribution 9. 

Proof. According to Subsection 4.1 of [47], we have N p satisfies Hy- 
pothesis 1. Then, Theorem 1 yields the desired result. □ 

2.3. Proof of Theorem 1. Motivated by the nonlinear filtering theory, 
the proof of Theorem 1 combines characteristics of (1) with classical tech- 
niques based on Girsanov's theorem. We start by recalling the result on the 
existence and uniqueness of the regular strong solution of (3) given by [47] 
(see also [45]). 

Definition 2. Let C satisfy Hypothesis 1. Suppose that T is either 
[0, +oo[ or [0, T] with T G [0, +oo[. We say that the stochastic process (^t(O)teT 
is a strong solution of class C of (3) on the interval T (for simplicity, C-strong 
solution) if: 

• (¥?t(£))ieT is an adapted process taking values in h with continuous sample 
paths. 

• For any t G T, E||v3 t (£)|| 2 < E||£|| 2 , ip t (£) G Dom(C) P-a.s. and 
sup se[0 ,t]E||G> s (0l| 2 <oo. 

• P-a.s., for all t G T, 

ft °° ft 

(8) <Pt(£)=S+ Go7rcM0)<fa + £/ L k on c (<P.(S))<W?. 

Jo k=1 Jo 
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Theorem 2. Let C satisfy Hypothesis 1. Suppose that £ is a ^Q-random 
variable taking values in rj such that £ £ Dom(C) a.s. and K\\£\\q < oo. 
Assume that T is either [0, +oo[ or [0,T] whenever T £ [0,+oo[. Then, there 
exists a unique C -strong solution (tpt(£))teT °f (3)- In addition, for all t € T, 
we have the following: 

(i) E|MOI| 2 = E||£|| 2 . 
(ii) E||C^(C)|| 2 <exp(at)(E||Ce|| 2 + at(E||C|| 2 +/?)). 

We now point out some immediate consequences of (Hl.l), (HI. 2) and 
(HI. 4) (see [47] for more details). 

Remark 2. Let assumptions (Hl.l) and (HI. 2) hold. Using the closed 
graph theorem, we see that G can be considered as a bounded operator 
from (Dom(C), (■, -) c ) into fj. By (2), for any x E Dom(C), £fcLi \\L k x\\ 2 < 
K\\Cx\\ 2 c , with K>0. 

Remark 3. Suppose that conditions (HI. 2) and (HI. 4) hold. Then, for 
any x in Dom(C), liuin^^ CP n x = Cx. It follows that 

Dom{C) = {x £ f) : (CP n x) ne z + is a convergent sequence}, 

because C is a closed operator. 

Remark 4. Let C satisfy Hypothesis 1. Applying Remarks 2 and 3, 
we see that G o ttq and L& o 7r<7, with k 6 N, are *8(Fj)/93(f))-measurable 
functions. 

The following proposition constructs a C-solution of (1) on [0, T] whenever 
T £ [0,+oo[ with the help of (3). 

Proposition 1. Suppose that hypotheses of Theorem 1 hold. Let (<Pt(£))t>o 
be the C-strong solution of (3), where £ is distributed according to 6. Define 
Q= \\ip T (£)\\ 2 -F, where Te]0,+oo[. For any t £ [0,T], we set 

x J>*(£)/IM£)II ; VMO^o, 
1 lo, tfMO = o, 

and 

(9) B? = W t k - f —±—d[W k M0]s, 

Jo WvsiOW 

with keN. Then (tt,3,($t)t€[o,T],Q,( x t)te{o,T],(Bf)^ T] ) is a C-solution 
of (1) with initial distribution 6. 
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Proof. Applying Ito's formula (or Lemma 2.1 of [47]), we obtain 

oo rt 



(10) IM£)I| 2 = II£H 2 + E/ 2Re(ip s (0,L Ws (0)dW t 



k 

s 




Consider the stopping time T n = inf{£ > 0: || </?*(£) II > n} AT, with n £ N. 
From Remark 2 and Theorem 2 we see that 

oo „T 

E E / Re 2 (<p,{Z),L k <p s (£))d8<n 2 K T (pU\\c + l), 
k =i Jo 

where Kt is a constant depending on T. It follows that 

' °° rt/\T n 



Y,] q 2Re(<p s (0,L k M0)dW* 



s) 

fc=i" u / te[o,T] 



is a square integrable martingale. Conditional Fatou's lemma now shows 
that (||^t(£)l| 2 )ie[o,T] is a supermartingale. Since E||(£ t (£)|| 2 = E||£|| 2 for all 
t>0, (|M£)l| 2 )te[o,T] is a martingale on (fi,ff, (fo) t6 [ ,T], F )- 

Let Q be the probability measure on J given by Q = ||</3t(£)II 2 ' ^i that 
is, Q is absolutely continuous with respect to P and the Radon-Nikodym 
derivative of Q with respect to P is ||¥?t(£)II 2 - The Girsanov-Meyer theo- 
rem shows that B k given by (9) is a continuous (O, #, (5t)ie[o,T]jQH oca l 
martingale. Since 

B k = W k - f t 2Re (<Ps(Z),L k <p s (Q} ^ 



' ' Jo [|^(0ll 2 

'1, if* = j 



foralltG[0,T],[^,^] t - 

According to Levy's theorem, B 1 ,!! 2 ,..., are independent Brownian mo- 
tions on (fi,$,(fo)t e [o,T],Q)- 

For all t £ [0,T], Q(v? 4 (£) = 0) = 0. Hence, combining (10) with (11) yields 

oo . t 



[12) 

oo rt 



vt(0f = Uf + Yl / \\M0\\ 2 ^ 2 (Xs,L k x s )d s 

oo . t 

+ E/ ll^(e)H 2 2Re(X s ,L fc X s )^. 
,, Jo 



k=i J0 



Let M 4 = EfcLi /o2Re(X s ,L fe X s )d5^. From Remark 2, Theorem 2 and 

EujRe (X s ,L k X s ) < E<Q\\L k X s \\ 

= E v \\L k <p a (S)\\ 2 , 
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we have (Aft)tg[o,ri is a continuous square integrable martingale. By (12), 

\\M0\\ 2 = U\\ 2 + fwMOfdiM+iM^i]^. 

Jo 
Therefore, ||v?t(£)ll 2 = exp(M t + [M,M]i/2)||£|| 2 , which implies 

yt{t)\\- l = eM-M t /2-[M,M} t /A)U\\- 1 . 
Hence, 



(13) 



k=i j ° 

E / IMOW^MXs, L k x s )dBl 
k=i Jo 



k= 
Using Ito's formula, (8) and (13), we obtain (4) with W replaced by B. □ 

Remark 5. Let assumptions of Proposition 1 hold. By the previous 
proof, the process (\\ft(0\\ 2 )t>o is a martingale on (0,5, (fo^g^T^P). Then 
for every T > 0, (||9?t(£)ll 2 )te[o,T] is a uniformly integrable martingale as both 
Section 4 of [34] and Section 4 of [35] stated without proof. 

Notation 1. Let E be a normed space. In the sequel, C([0,T],E) stands 
for the space of all continuous functions from the interval [0, T] to E endowed 
with the uniform norm || • ||oo. Moreover, we define C([0,T],M°°) to be the 
Cartesian product space YWLi C([0, T],R) equipped with the metric 

00 
d((f k ) k& , (g k )kem) = E 2" fe min{l, \\f k - g^}. 
k=i 

Proposition 2 deals with the uniqueness of solutions of class C for (1) in 
the sense of joint probability law when T is a bounded interval. 

Proposition 2. Let hypotheses of Theorem 1 hold. Take (Xt, (-B f fe ) feeN )t e [o.T] 
as in Proposition 1. If (Q, (Xt)t£[o,T]i (•^t)te[o,Tl) *s a C-solution of (1) with 
initial distribution 9 , then (X t ,(B k ) ke ^) t£ t 0T j and (X t , (-B t fc ) fcGN ) ie r 0i T] have 
the same finite- dimensional distributions. 

Proof. We begin by recalling that 

00 

fc=l 
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Then, (B k ) k ^ T , an( j (5 t fc )k|N are ranc j om variables taking values in C([0, T], 
E°°). Since C([0,T],M°°) is a Polish space, classical arguments (see, e.g., the 
proof of Proposition IX. 1.4 of [53]) reduce our proof to the case Xq = x and 
Xq = x for x £ Dom(C) with \\x\\ = 1. 

To deal with the above situation, we consider 6 = C([0,T], f)) xC([0,T],f)) x 
C^C^TjjIR 00 ) equipped with the usual product topology. For any / = 
if 1 J 2 J 3 ) G ©, we set Z k {f) = f k , provided k = 1,2,3. By the Yamada- 
Watanabe construction (see, e.g., proof of Theorem IV. 1.1 of [36] or proof 
of Theorem IX. 1.7 of [53]), there exists a probability measure [i on 03(0) 
with the following properties: 

• fio (Z 1 , Z 3 )- 1 = Q o (X, B)- 1 and (j, o (Z 2 , Z 3 )- 1 = Q o (X, 5)" 1 . 

• (^)te[o,T] is a sequence of real valued independent Brownian motions on 
(©, 0, (<5t)te[o,T] )/•*)• Here (0,(5,//) is the completion of (0,<B(0),/i) and 
for each £e[0,'r], 

<8t = fl ( ff (/( s ) :se[0,(i + e)AT])U{AG6: /x(A) = 0}). 

£>0 

In the sequel, j is either 1 or 2. From Remark 2 we have 

I'T 



JT /" ^z 2 (Zi,L k Zi)ds<K j (l + ^\\CZ{f)ds, 
k=1 Jo Jo 

and so X^Li Jo ^R^C^s'j-^fc^a) <^ s < °°- Then (i?t)t e ro ! Ti described by 

oo , t 
H t = -Y, 2Re(ZlL k Zi)d(Z 3 ) k s 

is a continuous square integrable martingale. For each n G N, we define 
T^=inf[tG[0,T]: / ||CZ^|| 2 ds > n| AT. 

Therefore, TjJ is a stopping time and [H n ,H n ]t < Km(n + T). Accord- 
ing to Novikov's criterion (see, e.g., Proposition VIII. 1.15 of [53]), we have 

(exp(i7 t n — [H Tn , H Tn ] t /2)) te [ 0T ] is a uniformly integrable martingale. For 
any fceN, we choose 

(14) W t k ' J = (Z 3 ) k + / 2Re(Z s \ L k Zl) ds. 

Jo 

Applying Girsanov's theorem, we deduce that (W t '"OjgroTi is a sequence of 
real valued independent Brownian motions on (0,(J5, ((Si) tG [ 0i T])/ i "')) where 

fi j = exp(Hp - [H T -,H T "] T /2) ■ /x. 
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Because (Z\ , Z 3 ) tG r 0i T] satisfies (4) with (X,B) replaced by (Z- 7 , Z 3 ) (see, 
e.g., Theorems 8.3 and 8.6 of [48] or Exercise IV. 5. 16 of [53]), we see that 

rtAT-' / oo \ 

(&)? =x+ n i G Zi + lY. Re 2 (Zi,L k Zi)Zi ds 

J ° \ fc=l ) 

rtATi / °° \ 

-J n2MZ J s ,L k Zi)L k Zijds 
+ £/ L k {Zi)dW^. 



k=V 



For any t G [0,T], we set rj t = ex.p(—H t /2 + [H,H] t /4). By, for instance, 
Theorem 11.37 of [50], 



00 rtKTi 



mj rv/\l n 

%" = 1-^ ri s Re 2 (Zi,L k Zi) ds 

k=i Jo 

oo ptAT-' 

+ W n VsRe(ZlL k Zi)dW^. 



The ltd formula leads to 

(Z^)f» = rr + / G(Z^)J» ds + £ / L k (Z> V )^ dW^ . 
Jo k=1 Jo 

Let (<fi( x ))te[o,T] be the solution of (8) described by Theorem 2 when 
the underlying filtered probability space is (0, &, { < St)te[o,T]i ^ 3 )- I n order 
to prove that (Z J ry tAT j)t e [o,r] and (ip 3 j(x))tg[o,T] are indistinguishable, we 
use the optional stopping theorem to obtain 

E^eMH? ~ [H T -,H T "] T /2)\& sATi ) 

where s £ [0,T]. Hence, 

r-tAT/ 



E^y o nc(z^n 2 d S 

=JfMW')^ll^ 



2 11,07.? ||2 

x E /X (exp(^" - [F T »,if T »] T /2)|© sATJ ))d S . 
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Then 

rttfi ■ ft 

E^J Q \\C(Z^)^fd S = j o E^l [0 ^(s)\\CZl AT ^)ds 

ptATl 

= E ll j Q l [0tTi] (s)\\CZi\\ 2 d S . 

This implies E^ l3 J tAT " \\C(Z j r))^\\ 2 ds < n. Similarly, E^||(Z^)^|| 2 = 1. 
Then, Ito's formula shows that fi J -a.s., for all t 6 [0, T], {Z 3 r\) t t j = if 1 3 (x). 

Since E^ J \\CZ ] t || 2 dt < oo, we see that 

(15) fi(H TJ - [H, H) TJ /2 = +oo) = 0. 
The integration by parts formula yields 

eJJ n 2Re(Zi,L k Zi)d(Z 3 )^j < oo. 

By (15), n is absolutely continuous with respect to fj? . It follows that /x-a.s., 
for all te [0,T], 

(16) ^)^Ti = < T ^)- 

Using again E M J \\CZ 3 t \\ 2 dt < oo, we have /i-a.s., for all t G [0, T], ^ > 0. 
Due to (16), we get fi-a.s. 

(17) [ t \\Crt(x)/yi(x)\\\\ 2 ds= f\\CZifds, 
Jo Jo 

for any t <T^. By [i(J ||CZ||| 2 ds = oo) = 0, 1 1— > J ||CZ^|| 2 (is is a continu- 

ous function //-a.s. Hence //(T^ = Tor J n ||CZ^ || ds = n) = 1. Applying (17) 
gives 

rT "||C^(x)/||^(x)|||| 2 ci S = n, 

provided T£ < T. In addition, (17) shows that j ||C^(x)/||^(x)|| || 2 ds < n, 
for any t <T£. Therefore, //-a.s. 

(18) Ti=mi^€[0,T\:J*\\C^(x)/\\ ( 4(x)\\\\ 2 ds>n^AT. 

Combining the Yamada-Watanabe construction with the uniqueness of 
the pathwise solution of (8) given in Theorem 2, we obtain 

M 1 o (^(x), W 1 )- 1 = n 2 o (<p 2 (x), W' 2 )- 1 
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on <B(C([0,T],f)) x C([0,T],R°°)). One way to see this is reasoning like in 
the proof of Theorem IX. 1.7 of [53]. From (18) we get 

{\\^{x)\\ 2 ■ n l )o^\ x ),W>\T l n )- 1 
(19) 

= (||^ 2 (x)|| 2 ^ 2 )o(^(x),W 2 ,T r 2)- 1 



on <8(C([0,T],f)) x C([0,T],R°°) x [0,T]). Since 

K^)H 2 = (v) 2 

= eM-H T j + [H,H] T3 /2), 
we see that \\tp J j (x)\\ 2 • ^ = fi. Then 

According to (14), (16) and (19), we have 

fi o ((Z 1 )^ , Z 3 )" 1 = /x o ((z 2 ) T », z 3 )- 1 . 

Because E^(/ ||CZs|| 2 <is) < oo, T^ /' n ~*oo T fi-a.s. Letting n — ► +oo 
yields n o (Z 1 , X 3 )" 1 = n o (Z 2 , X 3 )" 1 on <B(C([0,T], fj) x C([0,T],R°°) x 
[0,T]).Then,Qo(X,B)- 1 = Qo(A > , J B)- 1 on5S(C([0,T],f})xC([0,r],M oo )x 
[0,21). D 

We now obtain easily Theorem 1 combining Propositions 1 and 2 together 
with Kolmogorov's theorem. 

Proof of Theorem 1. Proposition 2 leads to the uniqueness of the 
C-solution of (1) with initial distribution on [0, +oo[. 

By Proposition 1, there exists a solution (P n , (-^f)tg[o,n]) (^t' n )te[o,n]) °f 
class C of (1) with initial distribution 9 on [0,n] for any n E N. Here 
(P,B'' n ) takes values in C([0,n],f)) x C([0,n],R°°). Let (X,S) be the co- 
ordinate mapping process defined on Q = C([0,+oo[, f)) x C([0, +oo[, R°°). 
Choose 5° = 58(0) and $ = o-((X s , B a ) : s < t). For every A G «8(C([0, n], f)) X 
C([0,n],R°°)), we define 

Q n (^- 1 (^))=P n ((X™, J B")GA), 

where vr n :0 -> C([0,n],f)) x C([0,n],R°°) is given by vr n (cj) = (u t )te[o,n] ■ 
Therefore, Q n is a probability measure on (0,$°). Moreover, Proposition 2 
leads to (Qn)neN being a consistent family of probability measure. Since 
(0,5°) is a standard measurable space, there exists a probability measure 
Q° on (0,5"°) such that Q° restricted on ^n coincides with Q n for every 
n G N. Let (3,(0) be the completion of (S°,Q°). Let (fo)t>o be the usual 
augmentation of (5"t)t>o- Using Proposition 1 yields (Q, (X t ) t > , (B t ) t > ) is 
a C-solution of (1) with initial distribution 9 on [0, +oo[. To see this, we can 
use, for instance, Theorems 8.3 and 8.6 of [48]. □ 
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3. Regular invariant measures. As a step toward understanding the large 
time behavior of the open quantum systems, this section treats Problem (ii). 

3.1. Main result. We first apply standard arguments to obtain the Markov 
property of the solution of (1). 

Theorem 3. Assume that hypotheses of Theorem 1 hold. Let f : (f), 53(f))) - 
(R, 23(R)) be a bounded measurable function. Suppose that (Q, (X t )t>o, (Wt)t>o) 
is the C -solution of (1). Then for any s,t>0, 

(20) E(f(X s+t )\$ s ) = E(f(X s+t )\X s ). 
In addition, we have 

(21) E(f(X s+t )\d s )= f f(z)P t (X s ,dz), 

where 

(99) p(v A ,_iQx(XfeA), xeDom(C), 

(22) r t {x,A)-\ SjB ^ x( £Dom(C\ 

provided A £ 23(f)) and (Q x , (Xf)t>o, (B' t ' x ) t >o) is the C-solution of (1) with 
initial data x G Dom{C). 

Proof. The proof is deferred to the Appendix. □ 

We now adopt a set-up similar to that in Section 2 of [24], where it is 
discussed as the existence of quantum stationary states. 

Hypothesis 2. Suppose that D satisfies Hypothesis 1 with C substi- 
tuted by D. Assume in addition that: 

(H2.1) The set {x S fj : ||-Da;|| 2 + ||x|| 2 < 1} is compact in the topology of fj. 
(H2.2) There exists a vector x belonging to Dom(D) such that ||x|| = 1 and 



E\\DXf\\ 2 ds < K^t, 



where K± is a constant depending on x and (Q,(Xf)t>o,(B' t ' x ) t >o) 
is the D-solution of (1) with initial data x. 

The next theorem provides a general criterion for the existence of a regular 
invariant measure for (1). 

Theorem 4. Let D satisfy Hypothesis 2. Then there exists a probability 
measure T on the Borel a-algebra of\), that is, 93(f)), such that: 
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(i) T(Dom(D) n{i£l): ||x|| = 1}) = 1. 
(ii) Let Pt(x,A) be given by (22) with C replaced by D. Then 

(23) T(A)=Jp t (x,A)T(dx), 

for any t>0, x€t) and A £ 53(f)). 
(iii) J h \\Dz\\ 2 T(dz) < oo. 

We can check directly condition (H2.2) of Hypothesis 2, for instance, in 
certain systems formed by an arbitrary number of identical Fermi parti- 
cles (see, e.g., Section 4.2 of [47]). Nevertheless, this assumption has the 
disadvantage of involving the solution of (1). This motivates the following 
hypothesis on the operators G and L& that guarantees the fulfillment of 
condition (H2.2). 

Hypothesis 3. The pair (C,D) of self-adjoint positive operators in f) 
has the following properties: 

(H3.1) The operator C satisfies conditions (HI. 2) and (HI. 4) of Hypothesis 

1. 
(H3.2) Dom(C)cDom(D). 
(H3.3) There exist constants (3 £ [0, +oo[ and N G Z + satisfying 



2Re(Cx, CP n Gx) + ^ \\CP n L k x\\ 2 < -\\Dxf + 0(1 + 
fe=i 



Ixll 2 ) 



whenever n> N and x G f) n . 
(H3.4) The operator D satisfies condition (H2.1) of Hypothesis 2 and con- 
ditions (Hl.l), (HI. 3) and (HI. 4) of Hypothesis 1 with C replaced 
byD. 

Remark 6. In many situations we can find an operator C such that 
the pair (C, y/aC) obeys Hypothesis 3 for some a > 0. This is the case, 
for instance, of the Jaynes-Cummings model of quantum optics and the 
multimode Dicke laser models (see Section 5 of [24] and Remark 2.1 of [47] 
for details). 

Remark 7. In Hypothesis 3 the function x i — ► ||C:r|| 2 plays the role 
of the Lyapunov function. In fact, the relation between the infinitesimal 
generator for an ordinary stochastic differential equation and the standard 
Lyapunov function (see, e.g., [1, 22, 38, 42]) is replaced by condition (H3.3). 

The theorem below provides an intrinsic sufficient condition for the exis- 
tence of regular invariant probability measures. 
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Theorem 5. Suppose that (C,D) satisfies Hypothesis 3. Then there ex- 
ists a probability measure T on 53(f)) for which properties (i)-(iii) of Theorem 
4 hold. 

It follows from the next lemma that condition (H2.1) of Hypothesis 2 can 
be expressed as the following: 

(H2.1)' D has finite-dimensional spectral projections associated with bounded 
intervals. 

Then the setting of Theorem 5 coincides essentially with the framework of 
Section 4 of [24] , where it is treated as the existence of quantum stationary 
states. 

Lemma 1. Let D be a self-adjoint positive operator in f) such that D is 
unbounded. Then conditions (H2.1) and (H2.1)' are equivalent. 

Remark 8. In Section 4, we apply Theorem 5 to a general model of 
harmonic oscillators with one mode. Moreover, [24] verifies essentially Hy- 
pothesis 3 in both the Jaynes-Cummings model and a multimode Dickc 
laser model. The analysis used in these three examples suggest to us that 
our main criterion for the existence of invariant measures is easy to check 
in a wide class of physical applications. Nevertheless, it is interesting to see 
how this criterion applies to other concrete quantum mechanics models. 

3.2. Proofs. We begin by establishing that the C-solution of (1) depends 
continuously on the initial data in a distribution sense. Theorem 6 is of 
independent interest. 

Theorem 6. Suppose that C satisfies Hypothesis 1. Let (x n ) ne ^ be a 
sequence of vectors of Dom{C) of norm 1. Assume that there exists x G 
Dom(C) for which Wmn^ooWxn- x\\ =0. Let (Q n ,(X^ n ) t > Q ,(B' t ' n ) t > ), with 
n G N , be the C-solution of (1) with initial data x n (i.e., Xq u = x n a.s.). 
Then for all t G [0, +oo[, 

(24) Q n o (X^)- 1 -> n=>00 Q o (Xf)~\ 

provided (®,(Xt)t>o,(Bt)t>o) is the C-solution of (1) with initial data x. 
Here the symbol =>• denotes weak convergence of probability measures. 

Proof. From Propositions 1 and 2 we see that, for any t G [0, +oo[, 



X n \-1 _ f\\,„ r„ M|2 tb\ „ f,„ („ \l\\.„(„ Mh" 1 

t 



o(x^y L = (\\Mxn)\\ 2 ^)o(Mxn)/\\Mx 



and Qo (X?)- 1 = (||^(x)|| 2 • P) o {(pti^/WM^W)' 1 - H ^e (ft(0)t>o is the 
C-strong solution of (3). 
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Let / : f) — > M be a bounded continuous function. Then 

^f{ip t {x n ))y t {x n )\\ 2 -^f( n {x))\\^(x)\\ 2 = H l n + Hl 
where H l n = E P /(^(x n ))(||^(x n )|| 2 - ||^(x)|| 2 ) and 

^ = E P (/(^(x n ))-/(^(x)))||^(x)|| 2 . 
Since (3) is linear, 

(25) E P ||(/5 4 (x n ) -93 f (x)|| 2 =E P ||^ f (x n -x)|| 2 < ||x n -a;|| 2 . 
Using the Cauchy-Schwarz inequality yields 

Hi < sup|/(2r)|(E P ||^(x n ) - ^(x)|| 2 ) 1/2 (E P (||^(x n )|| + |M*)||) 2 ) 1/2 
zefj 

< >/2sup|/(z)|(Ep||^(x n ) - ^(z)|| 2 ) 1/2 (|K|| 2 + ||(x)|| 2 ) 1/2 . 

zGh 

It follows that H\ ->0asn^oo. 

By (25), Po(i/?((x„)) _1 converge weakly to Po [ip t {x))~ l as n — > oo. Hence, 
there exists a probability space (£2, S^P) where a f)-valued random variable 
a (resp. a n ) is defined, with distribution Po (^(x))" 1 [resp. Po (<£>t(x n )) -1 ], 
and such that a n converge a.s. to a (see, e.g., Theorem 11.7.2 of [20] or 
Theorem 3.1.8 of [22]). Since E~||a|| 2 = E P ||y? t (2;)|| 2 < ||:r|| 2 , Lebesgue's dom- 
inated convergence theorem leads to H 2 = E~(/(a n ) — /(a))||a|| 2 — > n ^oo 0. 
Therefore, 

(26) (|MZn)|| 2 • P) (MXn))- 1 =^oo (IM*)|| 2 " P) ° M*))' 1 ■ 

For any x G f) , set 

x/||x||, if x / 0, 



vr(x) - 

1 0, it x = 0. 

Since (|Mx)|| 2 ■ P) ° (¥>t(aO) _1 ({°}) = 0, (26) implies 

((||^(x n )|| 2 .P)o(^(x n ))- 1 )ovr- 1 

^ n ^ 00 ((||^(x)|| 2 -P)o(^(x))- 1 )o^ 1 
(see, e.g., Theorem 1.5.1 of [12]). This becomes 

ff(irT^)d(\\^n)\\ 2 -r) 

J \\\<ft{Xn)\\ J 

J \\m{xn)V 

whenever / : rj —> K is a bounded continuous function, which establishes (24). 

□ 
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Remark 9. Let u = Dom(C) n {x G f) : ||x|| = 1}. The map from t) to 
E given by » i-> E Q:c /(Xf )l jDom(c)nu (x) + /(x)(l - l jDom (c)nu(a ; )) is mea- 
surable as soon as / : f) — > R is bounded and continuous. This follows from 
Theorem 6 and Remark 3. Then, a functional form of the monotone class 
theorem (see, e.g., Theorem 1.21 of [18]) yields that, for any A G 53(f)), 
x \—t Pt(x, A) is a measurable function from fj to R. 

Proof of Theorem 4. Let x be as in condition (H2.2) of Hypothesis 
2. Set u = {i£f): ||x|| = 1}. For every A G 5B(u), we define 



1 f n 
T n (A) = - / Q(X s G^)ds, 
n Jo 



where (Q, (X t )t>o, (-Bt)t>o) is the D-solution of (1) with initial data x. Using 
approximations by simple functions, we can assert that 

(27) fg(x)T n (dx) = -[ n Eg(X s )ds, 

Ju nJo 

the function g : (u, 2$(u)) — > (R, 58 (R)) being positive and measurable. 

Set Si = {x G u: ||-Dx|| 2 < /}, whenever / > 0. Combining Chebyshev's in- 
equality with condition (H2.2) of Hypothesis 2, we see that 



1 l" n 

r n (u\5,) = - / Q(||z?x s || 2 >o^ 

"-Jo 



1 /" n 
<— / E\\DX s \\ 2 ds 
nl Jo 

< K^l. 

It follows that (r n ) n£ N is a tight family of probability measures, since Si is 
compact. Therefore, there exist a subsequence (r nfe )/ cG p ! } and a probability 
measure T on 58 (u) such that T Uk converge weakly to T as k — ► oo (see, e.g., 
Section 1.6 of [12] or Theorem 11.5.4 of [20]). 

Let /„ im :u^R be given by f n ,m{ z ) = min {sup fc < n \\DP k z\\ 2 ,m}, with 
n, m G N. Using (27), condition (HI. 4) of Hypothesis 1 with D = C and 
condition (H2.2) of Hypothesis 2 yields 

f 1 f Tlk 

/ /n,m(z)r nfe (cfc:) = — / Ef nm (X s )ds 
Jv. n k Jo 

1 f Uk ,, ,,9 

< — / E\\DXJ\ 2 ds 



rik Jo 
<K £ . 

The function / njm is bounded and continuous, and so 

fn,m{z)T nk {dz) > k ->oo \ fn,m(z)T(dz). 



u 
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Hence, 

(28) [ fn,m( z ) T (dz)<Kz. 

Ju 
Suppose that sup neN ||DP n z|| 2 < oo. Then for any y G Dom(D), 
\(z,Dy)\= lim \{P n z,Dy}\ 

n — >qo 

< \\y \\ sup \\DP n z\\ 2 . 

neN 

Hence, z G Dom(D*) = Dom(D). By Remark 3, we have sup n6N ||D-P n z|| 2 < 
oo if and only if z G Dom(D). Therefore, for all z G f), 

t i \ t- tt \ v, t( \ \ \\Dz\\ 2 , if z G Dom(D), 

/„,»(*) /™ /(*) where /(*) = | ^ » . f ^ ^ ^^ 

Applying the monotone convergence theorem and using (28), we get 
J u f(z)T(dz) < K&. Thus, r(u \ Dom(D)) = and / u \\Dz\\ 2 T(dz) < K & . 
For every probability measure [i on 53 (u), we choose 



Ptti(A)= [ P t (x,A)(i(dx), 

Ju 



provided A G 53 (u). Let g : (u, 03 (u)) —>■ (R, 03(E)) be a bounded measurable 
function. Approximating 5 by simple functions gives 

(29) J g(x)P tf i(dx) = J U g(z)P t (x, dz)^ (x(dx) 

Since 

(30) T nk (un Dom(D)) = T(unDom(D)) = 1, 

the restriction of T Uk to 03(uPl Dom(D)) converge weakly to the restriction 
of T to 03(uPl Dom(D)). Here ufl Dom(D) is equipped with the relative 
topology induced on it by f) . Assume that g : u — > M. is a bounded continuous 
function. From Theorem 6 we see that the function from ufl Dom(D) to R 
given by x 1— ► J u g(z)P t (x,dz) is continuous. Then (29) becomes 

9 (z)P t r n , fe (dz) = / ( [ g(z)P t (x,dz))r nk (dx) 

unDom(D) JunDom(D) \Ju / 

-^k^oo ( g(z)P t (x,dz))T(dx) 

JunDom(D) \Ju / 

g(z)P t T(dz). 

unDom(D) 

From (30), we obtain that (see, e.g., Corollary 3.3.2 of [22]) 

(31) P t Tn k ^k-,00 PtT on u. 
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Using (21) and (27) yields 

1 r n k 
P t F nk (A) = — EP t (X s ,A)ds 
nk Jo 

= — E(E(l A (X s+t )\ds))ds, 
n k Jo 

where A £ «8(u). It follows that P t T nk (A) = ±£ k F(X s+t G A)ds. Hence, 
for any close set Finu, 

lim sup P t T nk (F) 

k-^oo 

= lim sup (r nk (F) + — ( P ¥(X S eF)ds- [ ¥(X S E F) ds] ) 

fc^ooV n k \J nk Jo J J 

= lim sup T nk {F). 

k— >oo 

Applying the Portmanteau theorem, we see that 

limsu P P t r nfc (F)<r(F). 

k— »oo 

Therefore, Pt^n k =^fc-^oo T on u. From (31) we have PtT = T. For any A E 
03(f)), we set T(A) = T(A n u), and the theorem follows. □ 

To prove Theorem 5, we need the following lemma which provides global 
estimates for E||Z?X t || 2 . 

Lemma 2. Let the hypotheses of Theorem 5 hold. Suppose that 9 is a 
probability measure on 55(f)) such that 9{Dom(C) D {x E 1) : \\x\\ = 1}) = 1 
and J^ \\Cx\\ 2 9{dx) < oo. Let (Q,(X t )t>o,(B t )t>o) be the C-solution of (1) 
with initial distribution 6. Then for all t > 0, 

(32) / E Q \\DX S \\ 2 ds < Eq||CX || 2 + 2(3t. 
Jo 

Proof. Assume that ^ is distributed according to 9. Let 4>-,n(0 be the 
continuous strong solution of the following stochastic differential equation 
on fj n : 

ft n ft 

(33) t/v(£) =P n i+ PnGljJs,n{® *S + V) / P„LfcV-,«(0 <*W* 

Applying Ito's formula yields (see, e.g., [47]) 

(34) E||<MOH 2 <E||£|| 2 . 
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Combining Ito's formula with condition (H3.3) of Hypothesis 3, we obtain 
that, for any n> N , 

<E\\cip , n m 2 +^[\-\msAO\\ 2 +{3a+UtAO\\ 2 ))ds. 

Jo 
Then, (34) becomes 

(35) ( E\\D^ n (Z)\\ 2 ds <E\\CP r £\\ 2 + 2(3T 
Jo 

for all n> N. 

Let v be the Lebesgue measure on *B([0,T]), where T is a positive real 

number. Because the unit ball of L 2 (fl x [0,T],P® v\ (Dom(D), (■, -)d)) is 

weak* compact, (34) and (35) imply that there exists a subsequence (ip., ni (£))/eN 

such that ip.,m(0 converge weakly in L 2 (£l x [0,T],P<8> u; (Dom(D), (■, •)£>)) 

as I ->• oo. From Section 3.2 of [47] we have E||i/>t )n ,(£) - ip t (£)\\ 2 -^/->oo 0, 

and so ip.,m(0 converge weakly to tp.(£) as I — ► oo in L 2 (Q x [0,T],P®!/; 

(Dom(D),(- r ) D )). Hence, 

/ E||^(0||l)^<hminf / E||^, ni (OI|?>d*. 

JO l^ooj 

By (34) and property (i) of Theorem 2, 

(36) [ E||^ t (0|| 2 dt<hminf / E\\D^ ni {0f dt. 

Jo l~*ooJo 

According to (35) and (36), we have ff M\\D(pt^)\\ 2 dt < E\\C£\\ 2 + 2/3T. 
Applying Propositions 1 and 2, we get (32). □ 

Proof of Theorem 5. Since Dom(C) C Dom(D), the C-solution of 
(1) with initial distribution 6 coincides with the .D-solution of (1) with initial 
distribution 8. Combining Theorem 4 with Lemma 2, we obtain the assertion 
of the theorem. □ 

Proof of Lemma 1. Suppose that assumption (H2.1)' holds. We define 
E to be the spectral decomposition of D (see, e.g., [55]). Then the following 
property holds. 

Property P. Every bounded interval /, with E(I) ^ 0, contains a point 
p such that 

(37) E({p})^0. □ 

In fact, there exists a nonzero vector y in the range of E(I) since E(I) ^ 0. 
We write E VtV for the finite positive measure (y,Ey). Take I\ = I. We define 
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the interval I n +i, with n G N, by the recurrence relation I n +\ = 1^"-, provided 
j n G {0, 1} satisfies 

(38) E ViV {lt)>E y>v {I^). 

Here I®, l\ is a partition of I n into two disjoint subintervals of the same 
length. 

Condition (H2.1)' implies that the dimension of the range of E(I) is equal 
to a natural number N. Then 1 = {n G N : Ey jy (I n ) ^ E y ^ y (I n+ i)} has at 
most N — 1 elements. Conversely, suppose that m, . . . , ra/v belong to X. By 

(39) Ey^y(l n ) = Ey 7 y(l n ) + Ey )y (l n ), 

for any n G N, we have Ey )y (I^~^ n ) ^ whenever n G {m, . . . ,nj^}. Since 

In?" 1 ,..., llT N , lt N N are disjoint, E(In7 Jni ), ..., E{ll~ ]nN ) and E(l'n N N ) 
are orthogonal to each other. It follows that the dimension of the range of 
E{I) is greater than or equal to JV + 1, which is impossible. Therefore, there 
exists iVo G N such that Ey t y(I n ) = Ey jy (lN ) for all n > No. 
Applying Lebesgue's dominated convergence theorem yields 

li n (t)E y Jdt) 




- Ey t y(l No ). 

Combining (38) with (39), we obtain E y>y (lN ) > 0. Thus, f] ne ^I n = {?}■ 
From 

\\E({p})yf = (y,E({p})y) 

= E y,y({p}), 

we have (37). 

Combining the definition of the resolution of the identity (see, e.g., Defini- 
tion 12.17 of [55]) with repeated application of Property P, we see that there 
exist real numbers pj, where j are natural numbers less than M S NU {+oo}, 
such that, for any k € N, the number of elements of {j :pj G [0, k]} is at most 
the dimension of the range of E([Q, k]) and E([0, k]) = J2pe\o,k] E({pj})- By 
the spectrum of D is concentrated on [0, +oo[, for any x G Dom{D) and 
y G rj , we have 

(y,Dx)= Y^Pj(yi E ({Pj}) x )- 

j<M 

Using the Banach-Steinhaus theorem yields Dx = J2j<AiPjE({Pj}) x (see, 
e.g., Theorem 12.6 of [55]). The operator D is unbounded, and so M = 
oo. Then p,- — > oo as j — > oo, which yields D satisfies condition (H2.1) of 
Hypothesis 2. 
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On the other hand, assume that condition (H2.1) holds. Then (D 2 + /) _1 
is a compact self-adjoint operator in I). Since is not an eigenvalue of 
(D 2 + /)~ 1 , there exist an orthonormal basis (e n ) ne ^ + of f) and a sequence of 
positive real numbers (/3 n )nez+ f° r which linin^oo [5 n = and (D 2 + I)" 1 = 
EnGZ + /?n(e n ,-)en (see, e.g., Theorem 19B of [62]). Hence, /3 n e]0,l], pro- 
vided n G Z + , and 



D = E V 1 //3n-l(en,-)en. 
ragZ+ 

This implies assumption (H2.1)'. D 

4. Application. This section illustrates the main results of this paper 
with a one-dimensional harmonic oscillator. 

Theorem 7. Let the assumptions of Example 1 hold. Suppose that p 
is a natural number greater than or equal to 4. Assume that 9 is a prob- 
ability measure on 53(f)) such that 9(Dom(N p ) n {x £ f) : ||x|| = 1}) = 1 and 
L ||./V p x|| 2 (9((ix) < co. Then, (1) has a unique solution of class N p with initial 
distribution 9 provided \a^\ > \a§\. 

Proof. Let (ej)j£Z+ be the canonical orthonormal basis on / 2 (Z+). 
Then Dom(NP) = {x £ t):J2]°=o f p \{ej,x)\ 2 < oo} and N p = E^oj^e,-, ')• 
It follows that N p is a self-adjoint positive operator in £ 2 (Z + ) and the con- 
ditions (HI. 2) and (H1.4) hold. Furthermore, the set {x £ l 2 (Z + ) : ||iV p x|| 2 + 
\\ x \\ 2 < 1} is compact since lim.j-+ 00 j p = oo. 

By Dom(N 4 ) C Dom(H) n flLiC 00 "*^*) n Dom(L* k L k )), the operator 
./V p satisfies hypothesis (Hl.l). 

For simplicity, assume that n > 2. A long easy computation shows that, 
for any x belonging to the linear span of eo, • ■ • , e n , 

6 n 

(40) 2Re(N p x, N p P n Gx) + ^ ||A^P n L fc 2;|| 2 < £cj-|sj| 2 , 

fe=i j=o 

where Xj = (ej,x) and 

cj = WiKMVl + VT+T)^- 1 + VjP2 P -2(j) + VJ+Tp 2p - 2 (j)) 

+ \ ai \ 2 (-2pj 2p +p{2p - l)j 2p ~ l + P 2p ~ 2 (j)) 
+ \a 2 \ 2 (2pj 2p +p(2p+l)j 2p ~ 1 +P 2p ^ 2 (j)) 
+ M 2 (-4pj 2p+1 + 8p 2 j 2p + P 2p -i(j)) 
+ \a 5 \ 2 (Apj 2p+1 + 8p(p + l)j 2p + P 2p -x(j)). 
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Here we use the same symbol Pi for different polynomials of degree I whose 
coefficients depend only on I. It follows that condition (HI. 3) holds when 
|« 4 | > I CK5 1 . Theorem 1 now leads to our claim. □ 

Theorem 8. Let the assumptions of Example 1 hold. Assume that p is 
a natural number greater than or equal to 4. Suppose that either \a^\ > \a^\ 
or |a 4 | = \a.§\ with \a 2 \ 2 — l a i| 2 + 4(2p + 1) |ct4| 2 < 0. Then there exists a 
probability measure V on 53(f)) satisfying (23) such that L \\N p \\ 2 T{dz) < oo 
and 

T(Dom(N p ) nj^Gt): ||x|| = 1}) = 1. 
Proof. We return to the proof of Theorem 7. Let |ck4| > \a§\. By (40), 

6 

2Re(N p x, N p P n Gx) + ^ \\N p P n L k xf 
fc=i 

n 

< K\x \ 2 + £ N V P (4pj(|a 5 | 2 - |a 4 | 2 ) + 2p (j)/j 2p ), 
i=i 

where K is a positive constant and n > 2. Here (0 2p (j))j e z + is a sequence 
for which lim.j^ 0O 02 P (j)/j 2p exists. For any r > 0, 

lim (rj(\a 5 \ 2 - |a 4 | 2 ) + 2p (j)/j 2p ) = -oo, 
and so the pair (C,y/~aC) obeys condition (H2.3) for any a e]0,4p(|a 4 | 2 — 

Kl 2 )[. 

Suppose that |a 4 | = \a$\ and \a 2 \ 2 — |ai| 2 + 4(2p+ 1)|«4| 2 < 0. From (40), 
we obtain 



2Re(N p x,N p P n Gx) + ^ \\N p P n L k xf 



k=l 



\ai\ 2 ) 



< K\x \ 2 + Y, \x,\ 2 j 2p {8p(2p + l)|a 4 | 2 + 2p(\a 2 \ 

+ 02 P (J)/J 2P ) 

where n > 2 and K > 0. Here (o 2p (j))j e z + is a sequence such that 

limo 2p (j)/j 2p = 0. 

Therefore, (C, \J~aC) satisfies the condition (H2.3) for any 
ae]0,2p(|a 1 | 2 -|a 2 | 2 -4(2p + l)|a 4 | 2 )[. 
From Theorem 5 we now obtain the claim of this theorem. □ 
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A particular case of Theorem 8 applies in the following simple damped 
harmonic oscillator. 

Example 3. In the setting of Example 1, consider 0\ = 0, fa = w, 03 = 
0, ai = v r My + 1), «2 = v3^ and ay = 0, (3 < j < 6). Here w, A and z/ are 
positive real numbers. 

Example 3 describes a mode of the quantized radiation field of an ideal 
resonator which interacts with two-level atoms that pass through the res- 
onator (see, e.g., [21]). In this situation, A is the energy decay rate, v is 
the number of thermal excitations in the steady state and uj is the natural 
(circular) frequency. Since \a 2 \ 2 — \oi\\ 2 = —A < 0, Theorem 8 gives the exis- 
tence of a regular stationary measure. This is in agreement with Section 4.1 
of [52], where the existence of a unique faithful stationary state is studied. 

Theorem 8 also covers the next basic radiation-matter interaction me- 
chanics. 

Example 4. In the context of Example 1, we define /?i, (3 2 , a±, a 2 , «3 
and a6 to be equal to 0. Moreover, we set /% € M, «4 > and as > 0. 

Example 4 simulates a two-photon absorption and emission process. Since 
the phenomenon of two-photon absorption was observed by Kaiser and Gar- 
ret in [37], models like Example 4 have been discussed in the physical lit- 
erature (see, e.g., [54]). Using Theorem 8 yields the existence of a regular 
invariant measure whenever 04 > as, which is in agreement with [27]. In [27], 
Fagnola and Quezada characterized all the invariant states corresponding to 
Example 4 with 04 > 05. 

APPENDIX 

We can prove the Markov property of the C-solution of (1) using tech- 
niques of well-posed martingale problems. 

Proof of Theorem 3. To prove (20), we modify the proof of Theo- 
rem 4.4.2(a) of [22]. Let s > 0. Consider the set A in &, such that Q(A) > 0. 
For any B £ J , we define Qi(B) = E(UE(1 B |&))/Q(A) and Q 2 (B) = 
E(1 A E(1 B \X S ))/Q(A). Then Qi,Q 2 < Q, that is, Qi, Q 2 are absolutely con- 
tinuous with respect to Q. 

Let r > 0. For any bounded measurable function g : (0, $ s +r) —* (M, *B(K.)), 
we have E((W s fc +t - Wg +r )g\$ s+r ) = 0, provided t>r and k G N. Hence, for 
any j = 1, 2 and fc£N, E Qj {{W* +t - W* +r )g) = 0, and so 

(41) E Qj (W* +t \$ s+r ) = W* +r . 
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Write B\ = W k +t - W t k whenever t > and k G N. Since Qj < Q, Q r a.s., 
for all t€ [0,+oo[, 

/s-\-t °° rs-\-t 

G(X r )dr + J2 L k {X r )dB k r . 

k=i Js 

Using (41), Qj <C Q and the Levy characterization of Brownian motion, we 
obtain that i? 1 ,!? 2 , . . . , is a sequence of independent Brownian motions on 
(ft,8 r ,(&+i)i>o,Q;)- Therefore, (0,5, ($ s +t)t>o,®j, (*.+t)i>o, (B k ) k g) is a 
(7-solution of (1) with initial data distributed according to the law of X s . 
By 

Qi o (Xs)' 1 = Q 2 o (Xs)- 1 = Q(X S G -L4), 

Theorem 1 leads to EQ 1 (/(A A s+i )) =E,Q 2 (f(X s +t)) for any t > 0. This gives 
(20). 

Let Q be the regular conditional distribution for (X, W) with range space 
C([0, +oo[, J>) x C([0, +oo[,R°°) given (X , W ). From Theorem 1 we deduce 
that Q-a.s. for all wESl, Q X()H o (X x °^\ B^M)" 1 = Q(w,-)- This fol- 
lows paraphrasing the proof of Proposition IX. 1.4 of [53]. Then for w in a 
set of probability 1 for Q, P t (X (uj),A) = Q(uj,tt~ 1 (A)), where A G 53(f)) 
and the map 7r is defined by ir(a,b) = a(t) for any a G C([0,+oo[, h) and 
6 G C([0, -r-oo^M 00 ). Using the definition of regular conditional distribution 
yields 

(42) E(l A (X t )\X ) = P t (X ,A), Q-a.s. 

for any AG 03(f)). 

Suppose that s is greater than 0. Set B? = W k +t - W k . Then (B^) k ^ is a 
sequence of independent (3s+t)t>crBrownian motions. Furthermore, we have 
(n,$,($ s+ t)t>o,®,(X s+t ) t >o,(B~£) k §®) is the C7-solution of (1) with initial 
data distributed according to the law of X s . According to (42), 

E(l A (X s+ t)\Xs)=Pt(Xs,A), Q-a.s. 

whenever i£!B(f)), and so (21) follows. □ 
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